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Abstract 

We derive the distributions of individual eigenvalues for the Hermitian Wilson Dirac Operator 
D5 as well as for real eigenvalues of the Wilson Dirac Operator Dw- The framework we provide 
is valid in the epsilon regime of chiral perturbation theory for any number of flavours Nt and 
for non-zero low energy constants We, 7,8- It is given as a perturbative expansion in terms of the 
/c-point spectral density correlation functions and integrals thereof, which in some cases reduces to 
a Fredholm Pfafhan. For the real eigenvalues of D\y at fixed chirality v this expansion truncates 
after at most v terms for small lattice spacing a. Explicit examples for the distribution of the 
first and second eigenvalue are given in the microscopic domain as a truncated expansion of the 
Fredholm Pfaffian for quenched -D5, where all fc-point densities are explicitly known from random 
matrix theory. For the real eigenvalues of quenched Dyy at small a we illustrate our method by the 
finite expansion of the corresponding Fredholm determinant of size v. 



1 Introduction 



At low energies Quantum Chrome-dynamics (QCD) has a very convenient analytical testing ground, 
the so-called epsilon regime of chiral perturbation theory (XPT) [TJ. In this limit standard XPT 
is reduced to a group integral over the zero-momentum modes that are then non-perturbative and 
dominate the theory, times a Gaussian free field theory of propagating Goldstone bosons. Although 
being an unphysical limit the analytic predictions in this regime are useful and can be tested against 
numerical results from lattice gauge theory. Of particular interest for the establishment of chiral 
symmetry breaking and the effect of gauge field topology is the Dirac operator spectrum. There 
are several equivalent methods to compute correlation functions of Dirac operator eigenvalues in this 
regime. Either one starts from XPT in a partially quenched also called graded setting which is based 
on supergroup integrals, see e.g. [2], or equivalently one can use replicated partition functions [3]. Yet 
another method that was actually first proposed in [4] uses chiral Random Matrix Theory (RMT), 
being equivalent to the epsilon regime of XPT at leading order (LO) [21 [5]. This latter approach is 
best suited to compute the distribution of individual Dirac operator eigenvalues [U [7J - which in turn 
are particularly clean objects for a comparison with lattice data. 

In the continuum limit this setup has already been very useful, in order to determine the low energy 
constants (LEC) in the chiral Lagrangian or to test lattice algorithms for their topology dependence, 
and we refer to [8] for a recent review including references. In XPT there are only two LEC to LO, 
the chiral condensate X which is the order parameter of chiral symmetry breaking, and the Pion 
decay constant F. Individual Dirac operator eigenvalues allow for a very detailed error analysis in 
the determination of S compared to the average density through the Banks-Casher relation, and we 
refer e.g. to [9] for recent work. In order to make the epsilon regime also sensitive to F a coupling to 
imaginary iso-spin chemical potential [i has been proposed and realised in [10]. This leads to realistic 
values of F when computing and optimising higher order corrections |llj . The computations were done 
in a partially quenched setting, where a /i-dependent Dirac operator is probed with /i-independent 
configurations. A wealth of results exists for the density correlation functions with [101 f!2l [TT] and 
without \x [H [13] , as well as for individual eigenvalue distributions with [TTJ [15] and without n [H [7J . 

What happens if one wants to include the effect of finite lattice spacing a, which is unavoidably 
present in all simulations? Perhaps surprisingly the effect is much simpler in Wilson XPT developed 
in |16[ [TT] , than in its staggered version [T8] . This is despite the fact that the Wilson term completely 
breaks chiral symmetry whereas the staggered Dirac operator partly preserves it. In a Symanzik type 
expansion Wilson XPT requires three extra LEC W6,7,8 to LO, that have to be determined in order 
to quantify the effect of a finite lattice spacing^, and we refer to jT9~] for review articles. 

The epsilon regime with Wilson fermions was put into focus in [20] , addressing mesonic correlation 
functions. Then two years ago the microscopic spectral density of the Wilson Dirac operator in 
the epsilon regime was computed [21] . The study of this operator has since seen a remarkably fast 
progress. First of all one has to distinguish between the non-Hermitian Wilson Dirac operator Dyy and 
its Hermitian version D§, and their corresponding densities. While we refer to ref. |22] for a discussion 
of the relation between topology and chirality v at finite-a, and to ref. [23] for the distinction of left 
and right eigenvalues for Dw, let us give a short list of what is known to date about the eigenvalue 
correlation functions. 

The quenched microscopic density for D§ and for the real eigenvalues of Dy/ were computed in 
|21[ [22] , respectively. These results were extended to Nf = 1 flavours in |24[ [25] and to Nf = 2 in 
[26j. There, also the quenched two-point eigenvalue density was given explicitly, and a computational 
scheme for higher densities and Nf > 2 was set up although it rapidly becomes cumbersome. In [27J 

1 This is so even though they represent unphysical lattice artefacts, with their values depending on the lattice action 
considered. 
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all quenched higher density correlations were computed for D5 at v = 0, 1. The corresponding results 
for the complex density of D\y eigenvalues including all v were given in [23J, where further results to 
all orders were announced which have just appeared [28J. Most of the previous results for the density 
(except |25j) were given for W% 7^ and Wqj = 0, which is expected to be a good approximation 
for a large number of colours [29]. However, it was also understood in [22 1 how to derive Wqj 7^ 
results from Wq 7 = results on the level of XPT by one extra Gaussian integration per LEC. This 
was now extended to an RMT setting in [28]. Furthermore, in [22] it was shown that in the epsilon 
regime QCD inequalities lead to a restriction on the sign of Ws to be positive for Wej = 0. This 
constraint was very recently extended to different combinations of LEC |30[ [3T| [32] . We note that the 
signs (of combinations) of these LEC are not merely of academic interest, but decide crucially on the 
scenarios of possible phase transitions, see e.g. the discussion in [TTJ, [301 E2]- Finally we would also 
like to mention the computation of the mean spectral density of -D5 in |33j . 

Based on previous experience at a = how useful individual eigenvalue distributions are, in [34"ll35] 
first lattice simulations to compare to the above epsilon regime predictions at a 7^ were made. There 
the individual eigenvalue distributions were generated numerically from RMT. We fill this gap here 
by providing the framework of how to analytically compute these distributions. This setup holds for 
general Nf and non vanishing LEC Wqj^. Our explicit examples are all quenched and at 1^6,7 = 0, 
but our framework can be used for all higher density correlations for D5 and Dy/ at Nf > that are 
available. 

The content of this paper is organised as follows. In section [2] we set up the general framework 
of how to compute individual eigenvalue distributions or their cumulative distributions in terms of 
density correlation functions which we assume to be given. In section [3] we then specify the joint 
probability distributions for D5- (subsect. 13. 1 1) and real D^-eigenvalues (subsect. 13.21) based on RMT. 
At Wqj = this leads to a closed Fredholm Pfaffian expression for the former, and a distribution 
of finitely many v real eigenvalues at small spacing a for the latter. Section [4] is devoted to explicit 
examples in the microscopic or epsilon regime. For quenched D5 we give the first and second eigenvalue 
in a truncated expansion of the Fredholm Pfaffian (subsect. 14. ip . In the case of Dw (subsect. 14. 2\i the 
approximate distributions are given in a finite series of expansion terms for v = 2,3. Our conclusions 
are summarised in section [5j 

2 Individual eigenvalue distributions: general theory 

In this section we describe the general formalism how to compute the distribution of individual eigen- 
values and their cumulative distribution, having the Wilson Dirac operator and its Hermitian coun- 
terpart in mind. We will focus on real eigenvalues here, although this setup can also be generalised to 
complex eigenvalues after choosing a family of contours to order them, and we refer to |36] for details 
in the case of a real chemical potential. 

Our starting point is chosen to be quite broad, given in terms of a partition function and its joint 
probability distribution (jpdf) of all eigenvalues that is invariant under permutations of eigenvalues, 
without further symmetry requirements. In particular we do not require the jpdf to be a determinant 
or Pfaffian in this section. As an example the jpdf can be obtained from the corresponding random 
two-matrix models and their representation in terms of D$- or real .Dyi/ -eigenvalues, which we will 
specify in section [3l Note that when introducing Wqj ^ by extra Gaussian integrations [22j [32j [28] 
the resulting jpdf is no longer a determinant or Pfaffian. Therefore it is important to be more general 
here. Our final expressions for individual eigenvalues will only contain /c-point spectral densities and 
their integrals. For that reason this setup also applies to partially quenched XPT where these densities 
can be generated by introducing pairs of fermionic and bosonic source terms. The analogous statement 
for spacing a = was reported in [37] , 
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We begin by defining the partition function of the theory under consideration, 

/oo \ n n 
UdXi V jpdf (\i,...,X n ) ■ (2.1) 
-00/ i=1 

It is given by n real integrals over the so-called joint probability distribution function Vj p df, which is 
taken to be invariant under permutations of all eigenvalues, without further symmetry requirement. 
It may depend on external parameters such as quark masses, the number of flavours or the chirality. 
We also define the fc-point density correlation function, 

re! If f°°\ n - k 

R k (Xi, . . . , Xk) = J^ZT~fc)\~z, yj J d\k+i ■ ■ ■ dX n Pjpdf (Ai, . . . , A n ) , (2.2) 

normalised by the partition function. The simplest example is just the spectral density for k = 1. For 
k = n the density is simply the normalised jpdf itself, R n = Vj p df/Z. The Rk are supposed to be 
known in the following as all other quantities will be expressed in terms of them. 
The k-th. gap probability or cumulative distribution is defined as follows: 

E k (b,c) = ^ k y ^ (J ) dXL-.dxJj + J ) dX k+1 ...dX n V jpdf (X 1 ,...,X n ) , (2.3) 

for k = 0, 1, . . . , re. It is proportional to the probability of having exactly k out of re eigenvalues inside 
the interval [b, c], and n — k eigenvalues outside this interval. Note that we do not specify any ordering 
among the k or re — k eigenvalues respectively, nor do we specify how many of the n — k eigenvalues 
are less than b or larger than c (we assume b < c without loss of generality). The simplest example 
with k = gives the probability of having [b, c] empty of eigenvalues. 

A remark is in order here. In the simplest jpdf we have in mind here that describes the transition 
from the Gaussian Unitary Ensemble (GUE) to the chiral GUE in RMT, the eigenvalues live on 
the full real line R and cannot be restricted to the positive half line R+. In particular there is no 
chiral symmetry as for the chiral GUE, that for every eigenvalue Xj > there is also an eigenvalue 
at — Xj < 0. Only on average the distributions of individual eigenvalues may be symmetric with 
respect to the origin in some special case (y = 0). Apart from that the derivation here is only a mild 
generalisation of [37] which we follow closely. Note also that we will still count the eigenvalues from 
the origin onwards, e.g. by choosing [b, c] = [0, c] for the positive eigenvalues, as this is where the 
microscopic large-n limit will be taken eventually. 

Rewriting + J°° = J*^^ — J fe c , applying the binomial formula and using the permutation in- 
variance of the jpdf we can rewrite the k-th. gap probability as 

re' 1 / f c \ k 

«M = ^Wzii) dXl - dXk 

- k /_ i„\ / roo \ n-k-l / rc\ I 

x v 



t fc ) (|°°) n dx k+1 ...dx n r jpdf (x 1 ,...,Xn) 

n-k j , „ c x k+l 

= Zj( -1 ) l |T (/ 6 J dX x ...dX k+t R k+l (Xx,...,X k+l ) , (2.4) 

where for k = I = the first term in the last sum is unity. For the simplest example k = we thus 
have 

E k=0 (b,c) = dX 1 R 1 {X 1 ) + ^ J dX 1 dX 2 R 2 {X 1 ,X 2 )-^J^ dX 1 dX 2 dX 3 R 3 (X 1 ,X 2 ,X 3 ) 

+ ... , (2.5) 
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where we only show the first few terms in the sum containing n + 1 terms. Note however that when 
considering the number of real eigenvalues of Dw, we can set to a good approximation n = v at 
small a, and then we will only need a small and fixed number of terms. In particular for n = 3 
eq. (|2.5p would be exact, and for n = 2 we only have the first 3 termg^l For the Hermitian Wilson 
Dirac operator D§, n will become proportional to the number of eigenvalues and thus will be taken to 
infinity. Even in that case taking only the first few terms in the now infinite series eq. (|2.4p will yield 
an excellent approximation, as we will see in section 0] and as was known already for a = [37j . 
The fc-th gap probability can be most conveniently written in terms of a generating function, 

E(b, c; = 1 + Ys(-O l l d\x . . . d\i R^,..., A/) . (2.6) 

We then have 



E k (b,c) = (-l) k — ^ E(b,c;0 



, for k = 0,1,... ,n . (2.7) 



The probability p k (b,c) of having one eigenvalue at the upper boundary A = c, with k — 1 > 
eigenvalues in [b, c] and all other eigenvalues outside this interval is defined as 

Pk(b,c) = k (^j^(^J b ^j dAi...dA fe _i (J + J \ dX k+1 ...dX n 

x'PjpdfiXi, • • • , Xk-i, Xk = c, Xk+i, • • • , A n ) , (2.8) 

for k = 1,2, ... ,n. For example if we choose b = and k = 1 this gives the distribution of the first 
positive eigenvalue on M + . The Pk(b, c) can be obtained by differentiation from the corresponding gap 
probability as we will show below. Analogously we can define the probability qk{b,c) of having one 
eigenvalue at the lower boundary A = b, with k — 1 > eigenvalues in [b, c] and all other eigenvalues 
outside this interval: 

q k {b,c) = k (^j^(^J b ^j dX ± . . . dXk-i (^J + J ^ dX k+ i . . . dX n 

x'PjpdfiXi, ■ ■ ■ , Afc_i, Afc = b, Afc+i, . . . , A„) , (2.9) 

for k = 1,2, ... ,n. This definition is convenient when counting negative eigenvalues from the origin 
onwards. For example the first negative eigenvalue on M_ follows by differentiating Eo(b,c) with 
respect to b and by choosing b < = c. 

The explicit relation between individual eigenvalue distributions and gap probabilities are as fol- 
lows, where we restrict ourselves to the pk(b,c). The relations for the qk(b,c) are similar. We obtain 
from Ek(b,c) by differentiation 

d c E k {b,c) = k\ { Pk {b,c) - Pk+ i(b,c)) , (2.10) 

(with po(b,c) = 0) or, after taking the sum on both sides, 

1-1 

pi(b,c) = -^2 T} d c E k (b,c) . (2.11) 



fc! 

k=0 



2 For n = 1 with a single eigenvalue no new information is gained from Eo compared to the density itself, see also eq. 
(|2.12|) for the first eigenvalue. 
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In the case of the simplest example k = 1 and 6 = the expansion eq. (12.5P maps to 



pi(0,c) = -— E o (0,c) = i?i(c)- J d\R 2 (X,c) + ... , 

for the distribution of the first positive eigenvalue on M + after choosing 6 = 0, and 

p 2 (0,c) = l dXR 2 (X,c) + ... . 
Jo 



(2.12) 



(2.13) 



for the second positive eigenvalue. Note again that if we had n = 2 here, the expressions in eqs. (|2. 12[) 
and (|2.13p would be exact, without further corrections. 



3 Application to the Wilson Dirac operators D$ and Dyy 

In this section we specify to the two RMT for the Wilson Dirac operator D\y and its Hermitian 
counterpart -D5 for an arbitrary number of flavours Nf, and give their respective jpdf. In the limit of 
large matrices specified in the next section [5] they describe the epsilon regime of Wilson XPT. 



3.1 The jpdf and quenched /c-point densities of D 5 

We begin with the RMT for the Hermitian Wilson Dirac operator with Wqj = in the form as it was 
introduced in [27], apart from putting an explicit n-dependence into the weight function here. The 
partition function is defined as 



Nf 

Z 5 (m;z) = / dH J dW det[£> 5 + z f l N ] exp 



D R 



ml n W 
-ml nH 



n 



(1 - a 2 ) 



n 

2a 2 



TrH 2 



+ H 



(3.1) 
(3.2) 



where we have introduced the Hermitian random matri^l H = of size N x N with N = 2n + u 
and the complex random matrix W 7^ of size n x (n + u). Here v G N will be kept fixed and of 
order one, it plays the role of chirality [22]. The parameter a 2 G [0, 1] incorporates the influence of 
the lattice spacing times the LEC, af at W%, see eq. (|4.2j) for the precise mapping between RMT and 
XPT quantities. It also allows to interpolate between the chiral GUE (a = 0) and the GUE (a = 1) 
in RMT when setting Nf = m = 0. The real parameters m and zj denote the standard quark mass 
and sources for ip'ysip, respectively. The jpdf for the eigenvalues Xj of -D5 was derived in [27] to where 
we refer for details, and we only give the result valid up to an overall (a-dependent) constant: 



V 5 ({x};m;z) = 




nx~- 
2^2 



N f 



n + z f) ) exp 



-nm 2 (-2n + iV(l - a 2 )) 
2a 2 (l - a 2 ) 



Z 5 (m;z) 



X &-2n+v{{x}) Pfl<^,j<A ^ ;l<g<^ 



f 2n+v \ 

] { / dxj V 5 ({x};m;z) 
s i=l J -°° J 



F(xj — Xi\m) x 



q-l 



e o- A 



-XJ e ^ 



0, 



(3.3) 



3 Note that in the original proposal [55] instead of a full matrix H only two Hermitian matrices were added on the 
diagonal blocks, similar to eq. p. 151) . 
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Here An({x}) = Tl1>j( x i~ x j) denotes the standard Vandermonde determinant. We have also defined 
the antisymmetric weight function 



F{x; m) 



exp 



nx 2 (l — a 2 ) 



erf 



4a 2 
XyJ n(l 



2a 



+ 



rri\/n 



+ erf 



XyJ n(l — a 2 ) 
2a 



rri\/n 



(3.4) 



appearing inside the Pfaffian. It could also be written in terms of the generalised incomplete error 
function erf(x,y) (c.f. |23|). Note that the Jacobian from the diagonalisation of D§ is not a simple 
Vandermonde determinant to some integer power, as it is in standard one-matrix models. The addi- 
tional Pfaffian determinant containing part of the weight function F(x; m) is a feature shared by the 
partition function for the non-Hermitian Wilson Dirac operator eq. (13. 14ft [23J, as well as by other 
non-Hermitian RMT, see |38j . 

The problem of computing individual eigenvalue distributions in terms of densities is now specified. 
All /c-point density correlation functions have been computed explicitly in the quenched approximation 
Nf = for v = 0, 1 [27J. For higher v > 1 for Nf = 0, 1 and Nf = 2 the spectral density has been 
computed explicitly in (22J [2U [25J [26] , respectively. There, the supersymmetric or graded eigenvalue 
method was used starting from the partially quenched chiral Lagrangian, but the results agree because 
of universality. Unfortunately from the density alone being the first term in eq. (|2.12p we do not 
gain any new information regarding the first eigenvalue. Let us emphasise that this lack of present 
information is not a restriction of our method. With the more general results for all fc-point density 
correlation functions for arbitrary Nf and v now at hand [28], these can be simply inserted into our 
setup without modifications. 

We continue to recall some of the results of [27] . For Nf = and v = 0, 1 all fc-point density 
correlation functions can be written as the Pfaffian of a 2 x 2 matrix valued kernel K n (x,y), 



Rk(xi, ...,x k )= Ph<i,j<k[Kn(xi, Xj)) , K n (xi,Xj) 



In(xi, Xj) 
— S n {Xj , X>i 



Sn^Xi, Xj) 



D n (Xi, Xj) 

For finite-n and v = its matrix elements read (we refer to |27] for details, including v = 1) 



(3.5) 



n -21L 

e 2a 2 



K =0 (x,y) 
i u n =0 (^y) 



'j-i 



S n~ ( X ,V) = [<p2j-2{x)R2j-l{y) - <p2j-l{x)R2j-2{y) 



, x 2 +y 2 



R 2 j-2(x)R 2 j-i(y) - R2j-i(x)R 2 j-2(yj) 



— (hj-2(x)(f)2j-i(y) - (p2j-i(x)(t>2j-2(y)) - F(x 



y;m) 



(3.6) 



with the corresponding skew-orthogonal polynomials and their integral transforms being given by 



R 2 j(x) 
R 2 j+i{x) 

(j)j(x) 



f. d-a 2 )^ r, se -w. 



j\ (l-a 2 yV2^ f°° 



{x + 2iasY 



m 



n- 



1 



(-1)JV5F 

O 2 

ny 

dy e~ 1 2a?F(x — y; m)Rj{y) 



ds e y/n{x + 2ias)Lj 



n- 



(x + 2iasY 



(I -a 2 ) 



(3.7) 



(3.8) 
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In particular the jpdf, eq. (|3.5j) for k = 2n + u, can be written as a single Pfaffian for both v = 0, 1: 



P5O1, . . . , X2n +J ,; m; z = 0) = Z 5 (m; z = 0) Pfi< ij < 2 „ +I ,[i^„(a; i , Xj)] , (3.9) 

which extends to v > 1 |28j. Inserting eq. (|3.5p into the generating function for all gap probabilities 
eq. (|2,6p it becomes a Fredholm Pfaffian [39J and can thus be written in a closed form 

n \ l c 

E(b,c;0 = l + ZX-O'yr / dsi...dxiPfi< i j<i[lirn(a?i J a;j)] 
i=i 76 



= JJ(1 - £Aj) = Pf[l - m . (3.10) 

i=i 

Here the Aj are the eigenvalues of the matrix integral equation A/(x) = dyK n (x,y)f(y). We will 
come back to the perturbative evaluation of this Fredholm Pfaffian after taking the microscopic large-n 
limit in section HJ 

Finally we would like to discuss the inclusion of af at WQj ~ a| 7 / 0, following [22]. There, it was 
first understood how to include these terms on the level of the XPT partition function, by introducing 
two extra Gaussian integrations and shifting the source terms for masses m and z, respectively. While 
the RMT partition function matches this only in the microscopic limilQ, it is tempting to assume that 
an RMT for ag 7 7^ can be constructed at finite-n, which was recently achieved in |28j . The integral 
transformed jpdf is non-trivial in general, and special care has to be taken regarding the normalisation 
when computing spectral densities and individual eigenvalue distributions of such RMT (see e.g. f l()' 
for a similar discussion). 

In particular it leads not only to a shift of the masses, but partly of the eigenvalues which can be 
viewed as extra auxiliary masses in the graded framework |22| (c.f. [32 \ 128]). If we define an RMT 
jpdf including ag 7 < [32] 

^ / poo \ 2 _ v§ _ vr 

V 5 ({x};m;z; a(i ;a 7 ) = — / dy 6 dy 7 e ^ ™WV 5 ({x + y 6 }; m - y e - z - y 7 ), (3.11) 

167r|a 6 a 7 | Vi-oo/ 

and a corresponding RMT partition function we obtain for the /c-point density 



2 2 
y 6 y 7 



(Ho) dy 6 dy 7 e ^ ^ 



Rk(x 1 ,...,x k ;m;z;a 6 ;a 7 ) = — — — — : — Z 5 (m - y 6 ; z - y 7 ; a 6 = 0; a 7 = 0) 

Z^{m; z; a§; a 7 )lbn\aQa 7 \ 

xR k (x 1 + y 7 , . . . ,x k + y 7 ;m - y e ; z - y 7 ;a 6 = 0; a 7 = 0) , (3.12) 

where we explicitly display all arguments here. We have checked that the resulting microscopic density 
k = 1 agrees with that following from the resolvent in [22j in various cases. Because eq. (|2.4p is a 
liner map this translates as follows to individual eigenvalues, 

2 2 

y 6 y 7 



(iZo) dy 6 dy 7 e^ ^ 



p k (c,b;m;z;a 6 ;a 7 ) = — — — — ; — Z 5 (m - y 6 ; z - y 7 ; a 6 = 0; a 7 = 0) 

Z^(m; z; a^; a7)lb7r|a6a7| 

><Pk(b + y 7 ,c + y 7 ; m-y 6 ;z- y 7 ; a 6 = 0; a 7 = 0) , (3.13) 
making our framework directly applicable in this setting too. 



4 The quenched partition function has to be m- and z-independent for that, unlike the convention in [27| , 
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3.2 The jpdf of D w 

We now turn to the RMT for the non-Hermitian Wilson Dirac operator following [22] , first at Wq 7 



0: 



■4 V 



Nr 



dA I dB I dWjJ det[D w + m f l N ] exp nTrWW^ - —Tr{A 2 + £ 2 )J , (3.14) 
/=i 



D 



w 



aA W 
-Wt aB 



(3.15) 



The Hermitian random matrices A = A^ and B = B^ are of sizes n x n and (n + z/) x (n + u), 
respectively. Here the parameter a > is not restricted and the masses rrif can in principle be taken 
non-degenerate. The Dirac matrix Dw can only be quasi-diagonalised and we refer to |23| for details 
where the jpdf was computed. Choosing v > without loss of generality, the In + v eigenvalues of 
D\y consist of n — I complex conjugate eigenvalue pairs, of I real right eigenvalues and of I + v real 
left eigenvalues. Here I can take any value < I < n, and we have to sum over all possible sectors. 
The partition function can thus be written as follows [23] for Nf = 0: 



AV 



2n+u 



dz 1R . . . dz nR dz 1L . . . dz n+uL V 



(N f ,u) i 



W 



ZlR, • • • , Z n R, ZlL, ■ ■ ■ , Z n +v h) 



(3.16) 



2n+v 



dz 1R . . . dz nR dz 1L . . . dz n+u L A 2n +v({z}) 



x det 

Ka<n, l<b<n+u, \<c<v 



9c(z aR , z bL )5W (z aR - z* bL ) + gu(x aR ,x bL )5^ (van)^ {ybi) 



1=0 



n—l 



dz\... dZn-i 



2l+u 



wx 



Zl, 



dx 1R . . . dx m dxiL . . . dxi +vL 

, Z n -l,X\ R . . . , X[ R , zl,... , Z*_ h Xn, • • • , Xi +u l) ■ 



(3.17) 



Here we defined two different jpdf, with a fixed number of / real right eigenvalues V^l > an d without 

fixing V^ f ' U \ We note in passing that the latter jpdf in eq. (|3,16p can also be written in terms of a 
single Pfaffian [28]. The one- and two-dimensional delta-functions in the variables z a = x a + iy a inside 
the determinant assure that we sum over all possible sectors. This sum is explicitly given in the last 
equation. The respective weight functions are defined as [23] 



9C (zi,z 2 ) 

gR(xi,x 2 ) 



rfi 


4 


* 

-4 


4vra 2 (l + a 2 ) 


\ Z 1 


- z 2 \ 


rfi 


4 


* 


47ra 2 (l +a?) 


\ Z l 


- z 2 \ 



exp 



- exp 



n(xi+x 2 ) 2 n(yi-y 2 f 



4a 2 4 
n(xi+x 2 ) 2 n(xi-x 2 ) 2 



(3.18) 



4o 2 



xerfc 



Vnjl + a 2 ] 



Fl ~ x 2\ 



2a 



9i(x) = 



n 



2vra 2 



exp 



nx 
2^ 
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(3.19) 
(3.20) 



The fact that we have to distinguish between left and right eigenvalues of Dw leads to more possibilities 
of defining density correlation functions, compared to eq. (|2.2p . When defining the density by inserting 
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a delta-function into the jpdf we encounter a density of complex, of left and of right real eigenvalues, 
respectively. For example following [23J the density of complex eigenvalues, R\c, an d the density of 
real right eigenvalues RmR is defined as 



Rm R (x)5^(y) + -R lc (z) = (jj dz 1R . . . dz n+vL V^ f ' v) {{z})8^ {z - z 1R ) , (3.21) 

and likewise for the density of real left eigenvalues Rmi ■ It is clear by comparing to eq. (|3.17p that 
in principle both densities get contributions from all sectors I > 0. Note that in [23J a particular 
combination of these densities was defined, the so-called density of chirality: 

Rr x ( x ) - R mL(x) - Rirr(x) . (3.22) 

Higher fc-point correlation functions have to be split as well. For example the 2-point function will have 
3 contributions: to find 2 complex eigenvalues, one complex and one real, and two real eigenvalues; 
and so on for higher k. 

We would like to point out that if one is only interested in the correlation functions of real eigen- 
values as we are here, at least approximately one can obtain the situation of having a jpdf of only 
finitely many v real left eigenvalues in the large-n limit. Hence for this setting expanding the Fred- 
holm determinant to obtain the first individual eigenvalue truncates at exactly v terms. In |23| the 
probability to obtain N a< j_d more real (right and left) eigenvalues in addition to v was computed as a 
function of the rescaled spacing a = a^/n/2 and of v, see figure 4 in [23J. For example for a value of 
a = 0.1 there are less than 2% additional real eigenvalues for v = (and even less for v > 0), that is 
2% contributions from sectors with I > 0. Restricting ourselves to the sector I = we can thus define 
an approximate jpdf for the real eigenvalues valid for small a: 

T : 'w,R V \ x ^L,---,x v l) = (J^j dz 1 ...dz n V^J^Q(z 1 ,...,z n ,zl,...,zl,xi L ,...,x l/L ) . (3.23) 



We can now follow eqs. (|2.2p and (|2.4p with n = v there, and we only need to insert the fc-point 
densities of real left eigenvalues for up to k = u, with v being small and finite. Note that in this 
approximation we have that 

R Rx (x) « R mL (x) > R mR (x) . (3.24) 

This is consistent when comparing to figures 5 and 6 in ref |23j for small values of d. 

Regarding the introduction of non-zero values for Wqj the same method discussed at the end of 
the previous subsection applies after adding a 75-mass term in eq. (|3. 15|) . and we refer to [32j[28] for 
more details on the transformed densities. 



4 Examples for individual real eigenvalues in the microscopic limit 

In this section we will give explicit examples for individual eigenvalue distributions. Because RMT is 
only equivalent to Wilson XPT in the epsilon regime after taking the microscopic large-n limit, we 
directly give results in this limit and refer to the literature for details. We present only quenched 
higher /c-point density correlation functions as an example, which we need in our expansion. However, 
as it should have become clear from the previous sections our method is not restricted to this case, it 
could immediately be applied to unquenched densities too. As a further simplification we choose to 
work here with W% = Wj = 0. Again this restriction is not a principle one, and we have shown how 
to lift it following [22j. To keep everything as simple and transparent as possible we choose to work 
with Ws 7^ only. 
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The examples for the quenched Hermitian Wilson Dirac operator are presented in subsection 
14, H and the examples for the distribution of real eigenvalues of the Wilson Dirac operator follow in 
subsection I4.2L 



4.1 Hermitian Wilson Dirac operator D 5 

The computation of the microscopic large-n limit was discussed in great detail in ref. [27] which we 
shall not repeat here, recalling only the results. The only difference to ref. [27] here is the extra factor 
of 2n in the exponent of the weight function in eq. (|3.1|) . We have chosen this convention here, to be 
achieved by a simple rescaling of all matrix elements, in order to have the same microscopic rescaling 
of the eigenvalues x, masses m and lattice spacing a as for Dyy in the conventions of ref. [23J. 
The microscopic quantities denoted by hat are defined as follows: 

2nx , m = 2nm, a = — \/2na . (4.1) 



x = 2nx , m = 2nm , a = ■ 

The map to the LEC in the chiral Lagrangian is provided by 

x = xVT, , rh = mVTi 



a 2 = af at VW 8 , 



(4.2) 



where V is the volume, S the chiral condensate, ai at the lattice spacing and W% the LEC. The 
microscopic fc-point density correlation functions have to be rescaled appropriately with n, 



Ps,k( x i 



,x k ) 



lim - — -r 

n-+oo (2n) ft 



Rk(xi, ...,x k ) 



(4.3) 



with all the arguments on the right hand side rescaled according to eq. (14. ip . The rescaled microscopic 
gap probabilities are defined as 

E s , k {b,c)= lim E k (b,c) , (4.4) 

n— ¥00 

and correspondingly the resulting rescaled individual eigenvalue distributions ps,k and qs,k follow. 
Hence the boundaries of the interval [b, c] considered have to be rescaled in the same way as the 
eigenvalues in eq. (14. lj) . 

We will now give the building blocks of the limiting microscopic kernel. All 3 kernels in eq. (|3.6p 
can be expressed through the kernel D n {x, y), the weight function F(x) eq. (|3.4p and integrals thereof, 
plus some extra terms for v = 1. Referring to |27] for details we have for our first example at v = 0: 

roo \ 2 



1 



32v / 2a 5 7r 



exp 



-.21 



16a 



;,2 



dsdr _„2_„2 
e 

7T 



(y — z + 4ia(s — r)) 



dtl ( y/t(fh 2 -(z + Aira) 2 )) I ^ yjt{m 2 - (x + 4isa) 2 )) 



dz exp 



(x 



^2^ 



32a 2 



erf 



— z + 2m n 
4v / 2d 2 " / 



+ erf 



2m \ 



4V2a2 



x,y) 



dzFs(y-z)S» s =°(x,z) - Fs(x-y) 



(4.5) 

s C-:,y), (4.6) 

(4.7) 



where 



F 5 (£) = exp 



32a 2 



erf 



+ 



m 



+ erf 



in 



is the rescaled, microscopic weight eq. 



dtI (XVt)I (YVt) 



Ay/20 2 2x/2d 2 J \4V2a 2 2V2a 2 

. One of the integrals can be computed analytically, 

2XI (Y)h(X) - 2YI (X)h(Y) 



(4.8) 







X 2 -Y 2 



(4.9) 
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However, this is not always advantageous when plotting the densities. Prom these building blocks we 
obtain the following expressions for the densities: 

Ps,i{x) = S s {x,x) , 
PSfl{x,y) = S s {x,x)S s (y,y) + Is(x,y)D s (x,y) - S s (x,y)S s (y,x) , (4.10) 

and for general k 



PS,k{xi, ■■■ ,Xk) = Pfl<M<fc 



I s (xi,Xj) S s (xi,Xj) 
-S s (xj,Xi) -D s (xi,Xj) 



(4.11) 



PS,1, PS,k=l,2,3 




Figure 1: Top plot: example for the approximated first, second and third positive D$ eigenvalue: 
ps,i{0,x) (left yellow curve), ps,2(0,x) (middle green curve) and ps, 3(0, x) (right pink curve), from 
eqs. (|4.12|) - (|4.14]) respectively, in the quenched theory with v = 0. The parameter values chosen are 
rh = 1 and a = 0.25. We also show the spectral density ps,i(x) (top blue curve) from eq. (|4.10|) and 
(|4.6j) for comparison, which is symmetric around the origin. Bottom plot: the spectral density eq. 
(|4.17[) and exact individual eigenvalue distributions (eq. (|4.15[) . [7J) for D5 at a = 0, shifted according 
to eq. ([4TT8]) . 

As our first example we will now give the distribution of the first, second and third positive 
eigenvalue of D5 for v = 0, using only the first 3 densities Ps~j° with j = 1,2,3 as an approximation 
when truncating the Fredholm Pfaffian expansion. In order to describe the first positive eigenvalue 
from the origin we choose [0, c] as the interval for the gap probability, differentiate with respect to the 
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upper bound c and rescale 2nc = x. Using eq. (12. lip at this level of truncation we obtain on 



Ps,i(0,x) 
P5,2(0,x) 
P5,3(0,x) 



rx -y rx rx 

PS,\{&)- J dyp s ,2(x,y) + - 



x rx rx 

dyps,2(x,y) - I 
o Jo Jo 

-y rx rx 

- / dydzp St3 (x,y,z) . 
1 Jo Jo 



o Jo 

dydzp Si3 (x,y,z) , 



dydzp s ,z{x,y,z) , 



(4.12) 
(4.13) 
(4.14) 



Their generating cumulative distributions are obvious. Note that when summing up these first 3 
approximate individual eigenvalue distributions eqs. (14. 12ft - (|4.14[> we trivially get the exact density. 
Would we add the exact individual eigenvalue distributions instead, this identity would only hold after 
summing up all eigenvalues, ps,i(x) = Yl'jLoPS,j( x )- 

Inserting the above expressions eqs. (|4.10p and (|4.1ip together with eqs. (|4.5p - (|4.8p leads to the 
curves shown in figure [T] (top plot). From it we see that our approximation gives well localised curves 
for the 1st and 2nd eigenvalue. While the 1st eigenvalue smoothly touches zero at around x ~ 6, 
clearly the approximation gets worse for the 2nd and even more so for the 3rd eigenvalue. To get a 
similar quality for the 2nd eigenvalue we would need to include up to terms containing ps,A, and for 
the 3rd eigenvalue terms up to ps t 5- The fact that the left tail of the first positive eigenvalue is not 
shown for this value of a is an artefact of counting eigenvalues from the origin. If we were to use 
another interval (say [— 2, x]) we could also determine this tail. However, the origin is a good point to 
choose since in the limit a — > chiral symmetry is restored. 

It is clearly seen when the approximation breaks down at latest: that is when p$ t i increases again 
from zero for large x, when ps t 2 becomes negative for the second eigenvalue, and when p$,3 exceeds the 
spectral density which is an upper bound to all individual eigenvalue distributions. This breakdown 
happens roughly at the same point x c ~ 7, up to which we can trust our approximated individual 
distributions (see also fig. [3] left). This is because from this point onwards the 4th eigenvalue would 
start contributing which is zero to our order of approximation. Nevertheless even the first part of the 
ascending curve of the 3rd eigenvalue could be used for fitting purposes to this order. In order to 
illustrate the described convergence of our approximation scheme we show in fig. [2] the same plot as 
fig- H (top plot), using one term less in the approximation: the density and the two-point function. 

Let us also compare our approximate distributions in fig. [I] (top plot) with the corresponding first, 
second and third eigenvalues when the lattice spacing a = 0, as shown in fig. [T] (bottom plot). In this 
limit also the individual eigenvalue distributions and microscopic density are known exactly: 



a=0 



a=0 



a=0 



2 6XP 



exp 



x 

'T 

x 2 

"T 



(J v {x) 2 - J u ^i(x)J u+ i(x)) . 



(4.15) 
(4.16) 
(4.17) 



Here we only give the first eigenvalue distributions following [6] . The expressions for higher eigenvalues 
can be found in [7]. In order to compare to the D5 eigenvalues which contain the mass, cf. eq. (|3.2p . 
we need to shift the density p~s,i eq. (I4.17P as in [22]. 



PS,l{X)\a=0 



m 



- 7 PS,i \\/ x 2 — rh 2 ) ®{\x\ — rh) , 

; ' V / a=0 



(4.18) 
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PS,l, PS,k=l,2 




Figure 2: For comparison we show the same plot as in the top fig. [ljfor the 1st and 2nd eigenvalue, using 
only an approximation up the second term in eqs. f)4. 12|) - f)4. 13|) with ps,\ and ps,2 at otherwise equal 
parameter values. It can be nicely seen how the quality of the approximation improves systematically 
from this plot to the previous plot, moving the breakdown point from x c ~ 4.5 here to x c ~ 7 there. 

and the individual eigenvalue distributions accordingly. These shifted distributions are plotted in fig. 
Q] (bottom plot). It can be clearly seen that the first eigenvalue is most sensitive to the influence of 
for which we have the best approximation. 
We would like to compare our new analytical predictions for individual D§ eigenvalues with RMT 
and lattice data, following very recent comparisons [6A\ [35] where numerically generated curves for 
RMT were used. As an illustration we use the RMT and lattice data from fig. 8 in ref. [33], shown 
in our figure [3J left and right, respectively. Our approximate analytical curves for the individual 
eigenvalues match very well the numerically generated RMT histograms from [33] in the left plot as 
they should. In the comparison to the lattice data - these are predictions from other data in [34J 
and not fits to these curves - one can see a slight systematic shift to the left. Our analytic curves 
clearly are a very sensitive tool to fit lattice data. Compared to fits of the spectral density also used 
in |34t [35] , individual eigenvalues have an unambiguous normalisation to unity. 



PS,k=l,2 PS,k=l,2 




Figure 3: Left plot: Numerically generated first and second RMT eigenvalue distributions (histograms) 
with 10 6 matrices of size n = 100 from [33] vs. our analytical curves. Right plot: comparison of 
lattice data (histograms) from [33] with our analytical RMT curves. The data set is for volume 20 4 , 
Pl w = 2.79, where the predicted values for £ = 216, rh = 3.5, a = 0.35 determined from another 
volume are used, and we refer to [33] for a more detailed discussion. 
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PS,li VS,k=l,2i QS,k=l,2 
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Figure 4: Top plot: the spectral density vs. the first two positive and negative individual eigenvalues 
PS,k=i,2(0, x) and gs,fc=i,2(£, 0) respectively, for m = 3 and a = 0.25 at v = 1. Negative values for 
the second eigenvalues are clipped. Bottom plot: the same curves at a = 0, with the spectral density 
and exact first and second eigenvalue given (eqs. (|4.17p . (|4.16p . |7J), respectively. The location of the 
delta-function due to the exact zero eigenvalues shifted to — m is indicated by the vertical line. 

Let us now turn to our second example at v = 1. Here the building blocks for the kernel read [2 7) : 

1 



D 



x,y) + 



16a7rv 2ir 



cxp 



x 2 + y 2 
16a 2 



dse s Iq I y rh 2 — (y + 4isd) 2 



dr e r 



rh + x + Aira 
y/m 2 — (x + 4ird) : 



Ji^\/™ 2 - (£ + 4ira) 2 ) - (4.19) 



dzF s (x-z)^- 1 (^,y) + 



47ra 



dse s Iq ( y to 2 — (y + 4isd) 2 

(4.20) 



1 



x(.Fs(x — z)e — Fs(y — z)e 



2 



dz exp 



z 2 + m 2 
16a 2 



dse s I V m 2 — (z + 4zsd) 2 



(4.21) 



We simply have to insert these into eq. ()4. to obtain the densities, and into eqs. (|4.12p - ()4.14p to 
obtain the individual eigenvalues to our order of approximation. The corresponding curves are shown 
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in figured] (top plot). Note that the spectral density is no longer symmetric with respect to the origin. 
This is due to the broadening by a > of the single zero eigenvalue {y = 1) which is shifted to the 
location of the mass, x = —rh. The first and second negative eigenvalue is determined separately 
using the quantities eq. (|2.9p . that is by differentiating the gap probability with respect to the lower 
bound of the gap [6,0]. When comparing to the same curves at a = in figured] (bottom plot), the 
first ("non-zero") negative eigenvalue there becomes the second negative eigenvalue for a > in the 
top plot, and receives considerable corrections as well. 

4.2 Real eigenvalues of the Wilson Dirac operator D w 

For fc-point correlation functions of left (or right) real eigenvalues of Dw much less was known com- 
pared to the previous subsection for D$. Explicit results included only the quenched k = 1 and 
(k = 2)-point function for arbitrary v. More general results for all quenched fc-point functions have 
been announced [23J and appeared now |28j. For k = 1 it is either given as the discontinuity of the 
corresponding resolvent [22j, or more explicitly as a double integral in [23] . For k = 2 its resolvent is 
given as a 4-fold integral over a 4 x 4 determinant times some rational function [26]. As expected this 
result can be factorised and brought into a Pfaffian form [28] . in analogy to [41] . Plotting the integral 
in the form of J2B] would still be a formidable task. It would give us a good approximation to the first 
real Dw eigenvalue and to the ascent of the second one. 

Rather than awaiting more tractable results announced in [23] that have appeared after finishing 
this paper [28] . let us illustrate our method in the small a-regime. Here it was shown [22] that 
the microscopic limit of the spectral density R^ x {x) becomes a v x v GUE distribution (just as the 
broadening of the zero-modes of D§). We therefore expect that the jpdf will be given by that of the 
GUE for v x v matrices. This is based on the analogy to -D5, where the factorisation of the jpdf into 
a GUE part for the broadened zero-modes and a chiral GUE part for the remaining eigenvalues was 
shown in [27]. We thus can write approximately for I = and v > 1 

VwM x il, ■ ■ ■ ,x uL ) ~ det [K u (x iL ,x jL )] , (4.22) 
where the GUE kernel inside the determinant is given by 



K v {x,y) = 




(4.23) 



In this regime all microscopic eigenvalues are measured in units of 4d. The spectral density 

pmUx) = K u (x,x) , (4.24) 

is then simply the Wigner semi-circle for finite n = v. Higher /c-point density correlation functions 
follow easily from eq. (|4.22p . 

PmUxiL, . . . ,x kL ) = det [K u (x iL ,x jL )] , (4.25) 

l<i,j<k 

for k = 1, . . . , v. Based on these expressions, we can now give all individual eigenvalue distributions 
by a finite number of terms for the gap probability. In this case it is a Fredholm determinant of size 
v x v, which after expansion gives rise to v + 1 terms. 

In figure [5] we show the spectral density and the individual eigenvalues for v = 2, 3. As the gap 
interval we have to chose here (—00, c]. For v = 2 for example the distribution of the first and second 
eigenvalues given in eqs. (|2.12p and (|2.13p are exact, without the dots. For v = 3 we can use eqs. 
(14.121) - (|4.14p for the first three eigenvalues which are now exact, and so on for higher v. 
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PlRL, PS,k=l,2 PlRL, PS,k=l,2,3 




-4 -2 2 4-4 -2 2 



Figure 5: The approximated density of real D\y eigenvalues for v = 2 (left plot) and v = 3 (right 
plot). The spectral density is given by the finite-^ GUE semi-circle, with the individual eigenvalues 
obtained as described in the text. All quantities are measured in units 4a. . 



We mention in passing that for a jpdf with such a small and fixed number of eigenvalues there 
exists an alternative computation instead of the Fredholm determinant. The individual eigenvalue 
distributions could also be computed as integrals over the ordered set of eigenvalues, 



Pi(xi) 

P2{x 2 ) 



.''2 



dx 



3 • • • 



■1-2 



X v -l 



dx?, 



•t'2 



dx v Vw»{x\, ...,x v ) 
dx v Vw,M.{xi, ... ,x u 



Xu-l 



x\<...<x v 



X\<...<X V 



(4.26) 
(4.27) 



etc., and then normalising them appropriately. 

We note that the computation of the individual eigenvalue distributions based on eq. (|4.22p also 
applies to the v broadened zero-eigenvalues of D$, due to the approximate density to R^ x (x) « Rs(x) 
for small a close to x = —rh [22J. One simply has to shift the argument of the Hermite polynomials 
in eq. (|4.23p as x — > x + rh. 



5 Conclusions 

We have set up a framework to analytically compute the distribution of individual real eigenvalues, 
both for the Wilson Dirac operator Dyy and its Hermitian counterpart D5 in the epsilon regime 
of Wilson chiral perturbation theory. These distributions are given in a perturbative expansion of k- 
point density correlation functions and integrals thereof. For small lattice spacing and for the real Dw 
eigenvalues this expansion truncates after v + \ terms where v plays the role of chirality on the lattice. 
As an example for quenched -D5 all fc-point density correlation functions are presented for v = 0, 1 
and we have given the resulting distributions of the first two positive and two negative eigenvalues as 
a function of the low energy constant W% times the lattice spacing squared as an illustration. Our 
results very well describe previous numerical results in the literature to the given order. We have 
briefly discussed how to include the effects of Wq^ in our framework. Hopefully this setup will become 
useful in the precise determination of all low energy constants to leading order. This is especially 
so because more analytic density correlation functions have just become available after finishing this 
paper [28] . 
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It would be very interesting to try to understand the mismatch discussed in [34] between lattice 
data and the analytical distribution of a single real eigenvalue of Dy/ (which is equal to the density 
Rr x as i n that case v = 1). However, considering possible higher order effects in chiral perturbation 
theory is beyond the scope of this article. 
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